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Abstract 

Within the next few years experiments at RHIC and the LHC will seek to create in the laboratory a 
quark-gluon plasma, the phase of matter through which the Universe passed very early in its life. It is 
believed that the plasma will survive long enough to reach thermal equilibrium. I give an introduction 
to the formalism of thermal field theory, the combination of statistical mechanics and quantum field 
theory needed to describe the plasma in thermal equilibrium, in a way that tries to keep close to the 
physics it describes. 


Introduction 

Thermal field theory is a combination of quantum field theory and statistical mechanics. This means 
that it is both difficult and interesting. The reason that we study it is that we want to describe the 
quark-gluon plasma, the phase that matter is believed to take above some critical temperature Tc- 
Lattice calculations suggest^^] that is about 100 MeV, or lO^^K. In the plasma phase the quarks 
and gluons are deconfined; they can move rather freely through the whole plasma. This is the phase to 
which the universe evolved soon after the big bang, and before the end of the century experiments at 
the new collider RHIC will try to re-create it in the laboratory, by making gold nuclei collide together 
head-on and dump their kinetic energy into a small volume. Similar experiments, at much higher 
energy, are planned later for the LHC at CERN. 

There is an obvious question: if a plasma is indeed produced, how will we know it? As yet there 
is no simple answer. There are estimates!^], necessarily based on very crude non-equilibrium theory, 
that suggest that the plasma will survive for a time long enough that it reaches thermal equilibrium 
before it eventually decays back into ordinary matter. So far, it is only equilibrium thermal field 
theory that is well formulated, and my lectures concentrate on this. For more information, the book 
by Le BellacI^J nowadays is the standard text, though the older book by Kapusta^ is still valuable, 
as is the classic Physics Report by Landsman and van Weert^^l. As I want my description to stay as 
close as possible to physics I will develop the theory using operators rather than path integrals, and 
mostly I will use the so-called real-time formalism. In thermal field theory it is very important to pay 
attention to the physics: some of the clever mathematical formalism that has been developed is very 
powerful for calculating certain quantitities, but it is all too easy to misuse it if one does not keep a 
close eye on the physics it is supposed to describe. 

Because in relativistic theory particles are continually being created and destroyed, it is appropriate 
to use the grand ensemble formalism of statistical mechanics. The grand ensemble consists of a very 
large number of copies of the system under study. The energy and other conserved quantum numbers, 
such as total lepton and baryon number, vary from system to system in the ensemble, but their 
totals over the whole ensemble are fixed, and therefore their average values per system are fixed. If 
the system is large enough, fluctuations about the average values are very small, and so the grand 
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ensemble formalism does effectively achieve the desired physical situation in which the conserved 
quantum numbers are fixed for each system. In practice we do not specify the ensemble in terms of 
their values; rather we use related quantities, the temperature and various chemical potentials. These 
arise as Lagrange multipliers when we maximise the appropriate quantity in order to calculate the 
most probable configuration of the ensemble. This mathematics leads to the grand partition function 

Z = (1) 


from which the macroscopic properties of the system in thermal equilibrium are calculated. Here (3 is 
the inverse temperature, (3 = l/ksT, and usually we use units in which Boltzmann’s constant ks = 1- 
The system’s Hamiltonian is H and is a conserved quantum number, such as baryon number, with 
^ the corresponding chemical potential. In the case of several conserved quantum mumbers, jiN is 
replaced with Yha k-a^a- 

The states |f) are a complete orthonormal set of physical states of the system. In scalar field theory 
all states are physical and so 

Z = tr (2) 

which is invariant under changes in the choice of orthonormal basis of states. In the case of gauge 
theories there are unphysical states, for example longitudinally-polarised photons or gluons, which 
must be excluded from the summation in (1). So then 


Z = tr P (3) 

where P is a projection operator that removes unphysical states. The presence of P can make things 
more complicated, and so to begin with I will consider scalar field theory, where it is not needed. 

All the macroscopic properties of the system in thermal equilibrium may be calculated from Z. In 
particular, for a system that is so large that its surface energy is negligible compared with its volume 
energy, the equation of state is 

PV = T log Z (4) 

Because we do not have sufficient knowledge of each system of the ensemble to specify which system 
state |f) it is in, we use a density matrix p to describe the system: in thermal equilibrium 

p = (5a) 

The “thermal average” of an observable corresponding to an operator Q, that is its average over all 
the systems of the ensemble, is 


<Q>=tiQp = Z-^ tr (56) 

Notice that, throughout, all operators and all states are familiar zero-temperature ones. The tem¬ 
perature enters only in the exponential, which defines how the zero-temperature states are combined 
together to define the statistical-mechnical ensemble used to calculate the thermal averages of the 
zero-temperature operators. 


Noninteracting scalar bosons 

For some systems of bosons there is no conserved quantum number N and therefore no corresponding 
Lagrange multiplier, that is there is no chemical potential p. For example, in the case of a heat bath 
of photons there is no constraint on their total number. Then in the scalar-field-theory case Z is just 
tr . 
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Consider first the case where interactions among the particles are so small that, once the system 
has thermalised, they can be neglected. This is a reasonable approximation for a weakly-interacting 
system, and is also a necessary preparation for setting up thermal perturbation theory, In the absence 
of interactions, the energies of the separate particles are good quantum numbers. To begin with, 
quantise the system in a hnite volume V, so that the single-boson energies are discrete. We can 
choose as the basis states |z) of the system those labelled by the single-particle occupation numbers 

Ur'- 

\n1n2 ns ...) 

and the eigenvalues of the noninteracting Hamiltonian Ho are 


nici -k n 2 e 2 + ^.363 -k ... 


So the noninteracting grand partition function is 

^-/3{niei+n2e2 + ...) 


{n.} 


HE- 


r Tlr 


, —/3n^er \ _ 


nrvi 


_ p-/3er 


(6a) 


and 

log Zo = - ^ log (1 - ) 

r 


In the continuum limit 



and so the noninteracting equation of state is 


H = — log Zo = -T 


d^k 


log(l 




( 66 ) 


( 7 ) 


( 8 ) 


where = \/k^ -|- If the bosons have non zero spin, there is an additional factor gg corresponding 
to the spin degeneracy of each single-boson state. 

In the continuum limit we usually work with fields 

/ d^k 1 

^ + h.c. (9a) 

with 

[a(k), a^ (k')] = (27r)^2fe°(5^^^ (k — k') (10a) 

In the discrete case, we usually dehne 

[Ur, ^1] drs (^^^) 

If we sum this over r, the result is 1. But if we apply V f d^k/(27r)^ to (10a), the result is rather 
2k^V. That is (10a) and (10b) have definitions of the operators a differing by a factor V2k^V. We 
correct for this by defining the field in the discrete case to be 

</.(x) = + h.c. (96) 
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The average number of particles in the single-particle level e* is 


(alas) = Zq ^ ^ (m n 2 ns ... |e ajas|ni n 2 n-a ...) 

{nr-} 

' E"-= /fe) 


(11a) 


where / is the Bose distribution 

It is interesting to calculate the average of the square of the number of particles in the single-particle 
level e^: 

((4a.)t= E"T“"'’‘- = (/(t.))T/(e,)(/(e.) + l) (lie) 

^ ■' Us 

The average of the square is not equal to the square of the average, because the number of particles 
in the level fluctuates from system to system in the ensemble. 

We can now calculate a noninteracting thermal Green’s function, defined as the thermal average 

{T<Pix)m)o = Zo" j;(i|e-^^«r<))(x)<))( 0 )|i) (12) 

i 


Observe first that there is an operator identity 

T<P{x)<P{0) = {0\T(j)ix)(j)i0)\0)+ ; 0(x)()>(O) : (13) 

where : : denotes the usual normal product. The first term contributes to (12) just the usual zero- 
temperature Feynman propagator. To evaluate the contribution from the second, use the discrete 

case (9b) and so obtain double sums Ylr s terms ajog, a^Os, aja|. When we take the necessary 
expectation values, only the first survives, and then only for r = s. So, using (11) we find that 

(: ^{x)m ■■)o = yJ2 + h.c. (14) 

r 

Going to the continuum limit, we may write 

/ fjiU 

(2vr)4 

DT{k) = — -^ -4 27r(5(/c^ - m^) n{k^) 

+ le 

where n{k^) = /(Ifc*^!). The second term is the contribution from the heat bath; it contains the 
(5-function because so far the heat-bath particles do not interact and so they are on shell. 


Perturbation theory 

Suppose now that we introduce an interaction and again calculate {T(p{x)4>{0)). 4>{x) is now the 
interacting Heisenberg-picture field: it is the familiar operator of zero-temperature field theory. We 
can develop a perturbation theory along the same lines as for zero-temperature scattering theory, by 
introducing an interaction picture that coincides with the Heisenberg picture at some time to: 

A(t) = ^-i{t-to)H 
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where Hqj is the free-held Hamiltonian in the interaction picture. Dehne as usual 

= U{ti,t2) 

so that A{t) = U{t,to). Because U is constructed from ordinary zero-temperature operators, as usual 
we have ^ 

U{ti,t 2 ) = Tex.p ^-i ( 17 ) 

We need U{ti,t2) for complex ti and t2, so that we integrate t along some contour C running from t2 
to ti in the complex plane and generalise the time ordering T to an ordering Tc along C: the operator 
whose argument is nearest to ti along the contour comes hrst. 

Now, from the dehnition of U, 

[/(to - i/3, to) 

so that 


Z ^tr e (p{x)(p{0) = Z ^tr e ^^°^[/(to — i/3, x*^) (/>/(x)[/(x‘^, 0) (/>/(0) [/(O, to) 

= ZoZ-\U{to - i/3, x°) cPiix) U{x^,0)<j)i{0) [/(O, to))^ 

where I have used U{ti,t2)U{t2-,t^) = [/(ti,t 3 ). 

Compare with what one has at zero temperature; 

(0|(/)(x)(/)(0)|0) = ( 0 |[/(oo,x'^)(/) 7 (x)[/(x'^, 0 )(/>/( 0 )[/( 0 , —oo)|0) 

In (18) we have a non-interacting thermal average instead of a vacuum expectation value, the hrst 
argument of the hrst U is (to — i/3) instead of oo, and the second argument of the last U is to instead 
of —oo. In fact there is a close similarity between thermal perturbation theory and the usual zero- 
temperature Feynman perturbation theory. One draws familiar-looking Feynman graphs, now called 
thermal graphs, and the main difference is that instead of the internal lines in the ordinary Feynman 
graphs representing vacuum expectation values, in the thermal graphs they are non-interacting thermal 
averages. 

In order to derive this, one needs to establish Wick’s theorem. It is a remarkable fact that indeed, for 
example, 

{Tc(t)i{xi)4>i{x2)(t)i{xi)(t)i{x4))o = '^{T4>i4>i)o {T4>i4>i)o (19a) 

where the sum is over the possible pairings of the fields. For almost every ensemble other than one 
in thermal equilibrium there would be correction terms to (19a). To understand this, go back to the 
discrete case. One of the relations needed to establish (19a) is 

(aja|ata^^) = (ajat)(aja.u) (aja„)(ajat) (196) 

This relation is straightforward except when the indices r,s,t,u are all equal. Then, from (11), the 
right-hand side is just 2(/(er))^, while with (10b) the left-hand side is 

((aja,.) (aja,.)) 

which we calculate using (11c). We get the same result as for the right-hand side because of the 
fluctuation term /(/ -I- 1) in (11c); for most other ensembles this term would be different and so 
Wick’s theorem would not be valid. 

One needs to choose a value for to- A common choice is to = 0, with the contour C for the t integrations 
running from 0 to —if3 along the imaginary axis. This is the imaginary-time formalism. Alternatively, 
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to ^ —oo, which with suitable contour choice gives the real-time formalism. Although under certain 
circumstances the imaginary-time formalism is a powerful calculational tool, the real-time formalism 
stays much closer to the physics and is more versatile. 


Real-time formalism 

One is interested in equilibrium properties of the plasma at finite times. Presumably these are inde¬ 
pendent of how it reached thermal equilibrium. So, as in familiar scattering theory, we are free to 
imagine that the interaction slowly switches off as we go into the remote past, and then when we take 
to ^ —oo the interaction-picture fields become the usual noninteracting in fields. The corresponding 
in states are direct products of non-interacting single-particle states. We need to choose how the 
contour C runs from —oo to (—00 — if3), and the choice that keeps the formalism in the most direct 
contact with the physics is the so-called Keldysh one: along the real axis from —00 to 00 , back to 
— 00 , then straight down to (—00 — i(3): 




For most applications (though not all!®]) it turns out that the vertical part of the contour may be 
omitted. Then we may write the propagator that corresponds to a line of a thermal graph as a 2 x 2 
matrix; 


(r(/)in (xi) (/)in (X2 ) ) 0 (0in (X2 ) </>in ) ) 0 
((/)i^(xi)(/)in(x2))o (T(/>in(xi)0in(x2))o 


When both xi and X 2 are on the —oo to oo part of C, the ordering Tc along the contour is ordinary 
time ordering T; this corresponds to the element Du of D. When both are on the oo to —oo part 
of C, Tc is anti-time-ordering T; this corresponds to D 22 - The off-diagonal elements correspond to 
Xi being on one part of C and X 2 on the other. There is translation invariance; the elements of D 
depend only on the difference between xi and X 2 - 


I have already shown how to calculate Du] the result is given in (15). The other elements of D may 
be calculated in the same way, and its Fourier transform is 


D(fc) 


_ i _ 

2Tr5~^{k‘^ — m^) 


27rd (k^ — m^) 

— i 

k'^ —m'^ — ie 


-|- 27r6{k‘^ 


m^) n{k^) 


1 

1 


1 

1 


(21a) 


It may also be written in the form 


D(A:) = M(A:°) D(fe) M(A:°) 


(216) 


with 


M{k°) 

t){k) 


=V^n(/cO) 


'eh0\k°\ 


ei/3|fc°l 


D{k) 0 
0 D*{k) 


D{k) 


i 

k2 _ jyip, _j_ 


(21c) 


For the case of a fermion field, there is a rather similar matrix propagator, but with the Fermi-Dirac 
distribution replacing the Bose distribution. 
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Matrix structure 

The elements of the matrix propagator (20) are not independent. For example, 

D 2 i{x) = ((/)in(x)0in(O))o = e“^^°‘"(/)in(x)(/)in(0) 

= Zf^^tr (j){x^ — i/3, x) 

= T)i 2 (x° - i/3,x) (22a) 

Here, I have used a general property of traces, that tr(ylH) = tr(i?yl), and the fact that //oin is the 
time-translation operator for the noninteracting field (pin- The Fourier transform of (22a) is 

Di2{k) = D2iik) (226) 

Also, from their definitions (20), one can see that Du and D 22 may be expressed in terms of D 12 and 
D 21 . It is this, together with (22b), which is responsible for the matrix structure (21b). 

Define a dressed thermal propagator matrix D'(xi,X 2 ) analogous to D(xi,X 2 ) in (20), but with the 
interacting Heisenberg field instead of (pin- For example, 

D[ 2 {xi,X 2 ) = {(pix 2 )(pixi)) (22c) 


Then, because H is the time-translation operator for (p, we can again derive 

D[,{k) = D',,{k) 

and so deduce that D' has a matrix structure similar to that of D: 


D'(A;) = M 


D'{k) 0 

0 D'*{k) 


M 


(22d) 


(23) 


where M(A:*^) is the same matrix as in (21c). Define the thermal self-energy matrix 11 by 


-in = D“^ - D' ^ 


Then it follows that 11 has the structure 


-iU = M 

If we then solve (24a) for D', we find 


-1 


-fn(A:,T) 

0 


D' = M ( 

0 


0 

Hn(A:,r)]^ 


0 

— i 

fc2-m2-n* 


M 


-1 


M 


(24a) 


(246) 


(25) 


So it is natural to interpret Re H as a temperature-dependent shift to the mass m^. H also has an 
imaginary part, so the propagation of the field through the heat bath decays with time. 


In scalar held theory, 


-in = 



+ 



+ • • • 


To calculate the contribution to ni 2 from the second term, for example, one needs 
1 1 2 2 



12 12 
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where each line ' 
opposite in sign. 


j represents Dij{k) and each vertex 2 is the same as the normal vertex 1, but 


Imaginary-time formalism 

The real-time formalism stays close to the physics, but has the calculational complication that the 
propagator is a matrix. In the imaginary-time formalism there is not this complication, though except 
for a few simple cases there is the need to perform an anlytic continuation from imaginary to real 
time at the end of the calculation. In the imaginary-time formalism the t integration runs along the 
imaginary axis, so ordinary time-ordering is replaced with ordering in imaginary time; 


D{xi,X 2 ) = 6<(—Im t)D 2 i{xi,X 2 ) -|- 0(lm t)Di 2 {xi,X 2 ) (26) 

where t = x^ — x®. Because both x^ and x% are integrated from 0 to ij3, we need D{xi,X 2 ) for values 
of Im t in the range —(3 to +P- In this finite interval it has a Fourier-series expansion; 

OO 

= ^ (27a) 

^ n= — oo 


where uOn = mrjp. However, the relation (22a) implies that D{t,x) = D{t + i/3, x), so that only 
even values of n contribute to the sum. In the case of fermions, the anticommutativity of the fields 
results in a minus sign appearing in the corresponding relation (22a), and so then only odd values of 
n contribute. 


If we apply a 3-dimensional Fourier transformation to (27a) and invert the Fourier summation over n, 
we find 

/•i/3 


DM= / dte-^-^D2i{t,]^) 


(276) 


JO 


which turns out to be just the ordinary Feynman propagator with = iiVn- So the Feynman rules 
are just like the zero-temperature ones, except that the energy-conserving (5-function at each vertex 
is replaced with a Konecker delta which imposes conservation of the discrete energy, and round each 
loop of a thermal graph 



(28) 


Gauge theories 

For gauge theories there is the complication that the grand partition function has to include the 
projection operator P onto physical states; see (3). The manipulations that set up the scalar-field 
perturbation theory include the use of the basic property of traces, trace AB =trace BA. To follow 
exactly the same route for a gauge theory one would need to use trace AB¥ = trace BAF, but this is 
not a valid identity. It turns outl^] that there are two formalisms for perturbation theory for a thermal 
gauge theory; 

1 Only the two physical degrees of freedom of the gauge field (the transverse polarisations) acquire 
the additional thermal propagator; the other components of the gauge held, and the ghosts, 
remain frozen at zero temperature. (This is for the bare propagators; self-energy insertions in the 
unphysical bare propagators do depend on the temperature.) 

2 All components of the gauge held, and the ghosts, become heated to temperature T. 

In the zero-temperature held theory, the ghosts are introduced in order to cancel unwanted contri¬ 
butions from the unphysical components of the gauge heldl®], and the two formalisms of the thermal 
perturbation theory lead to the same answers for calculations of physical quantities for that reason. 
Often, using formalism 1 makes calculations simpler. It also makes them stay closer to the physics. 



Photon or dilepton emission from a plasma 

As an application, consider the emission of a real or virtual photon of momentum q from a quark-gluon 
plasma. This is supposed to be an important diagnostic test of whether a plasma has been created in 
an experiment and has reached thermal equilibrium, and is then a way to measure its temperature. 

Before a photon is emitted, the plasma is described by the density matrix 

p = ^ |zm)(zin|e“^'^ (29) 

i 

The emission probability is calculated from squared matrix elements of the Heisenberg-picture elec¬ 
tromagnetic current!^]: 


= Z J (i"^xe*'^'^(/out| |im)(iin|e 7*^(0) |/out) (30a) 

where I have introduced also a complete set of out states for the plasma. These satisfy the completeness 
relation 

X] I/out) (/out I = 1 
/ 


and so 


Wf^’^iq) = J d^xe*«'^(iin|e-^^J‘^(0)J^(x)|zm) 


(306) 


If we introduce a matrix G^‘'{q) analagous to D', but with thermal averages of products of electro¬ 
magnetic currents instead of fields, IT^'^(g) is just G 12 ('i)- So we draw thermal graphs where current 
q enters at a 1 vertex and leaves at a 2 vertex, and distribute the labels 1 and 2 in all possible ways 
on the other vertices. For example, 



where the solid lines are quarks and the dashed lines are gluons. The emission rate is calculated 
from a sum of matrix elements times their complex conjugates; the vertices labelled 1 correspond 
to a contribution to the matrix element and those labelled 2 to its complex conjugate. The sets of 
1-vertices and of 2-vertices are joined by (12)-lines which, according to (21a) are on shell and represent 
particles in the heat bath. According to (21a), 

f) = i 27r6{k^-m^)n{k^) (k° > 0) . . 

( 27r(5(A:^ — m^) (1-|-n(/c‘^)) (k^ < 0) 

So, because the density of particles with momentum k in the heat bath is n{k^), when k^ > 0 a 
(12)-line represents a particle with momentum entering the photon-production reaction from the heat 
bath, while when < 0 it represents instead stimulated emission of a particle with momentum —k 
into the heat bath. In the loop integrations, the energy on each line extends over both positive and 
negative values, so the thermal graphs sum together many physical processes. Consider the first graph, 
for example. Its right-hand part, with the 2-vertex, represents the processes 





(ordinary Feynman graphs) 
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depending on the signs of the relevant internal energies. In fact, energy-momentum conservation allows 
only the first one to be non-zero. I have not drawn in the other heat-bath particles, but remember that 
they are there as spectators. The other part of the thermal graph contains 1-vertices and (ll)-lines, 
and according to (21a) the (ll)-propagator is the ordinary Feynman zero-temperature propagator, 
plus a thermal part. If I use only the zero-temperature part of Du in each of these lines, I obtain the 
amplitude 

(ordinary Feynman graph) 



(plus other terms which again vanish for kinematic reasons) and so part of the thermal graph represents 
the interference between this and . If instead I use the thermal part n{k^) 27r6{k‘^) of the (11) 
gluon propagator, I obtain the amplitudes 



(ordinary Feynman graphs) 


In each case, the incoming and outgoing gluon lines must have the same momentum k, so that these 
amplitudes again interefere with ^ but now with the gluon k being one of the spectator particles 
in the heat bath: 


k 


(ordinary Feynman graph) 



Similarly, I can identify physical processes that involve the thermal parts of the (11) quark propagators. 

Even a simple-looking thermal graph corresponds to a large number of physical processes^^'^l, each of 
which can be rather complicated. An example is 



for which just one of the physical processes is the interference between 


and 



(ordinary Feynman graphs) 


A disconnected Feynman graph occurs because the thermal graph has “islands” of groups of 1-vertices 
entirely surrounded by 2-vertices. 
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Resummation and discontinuities 

The matrix amplitude G^’^{q,T) has structure similar to that of D' in (23); 

(32) 

From the explicit form (21c) of M we find that 

Gi 2 ((?, T) = (g((?, T) + G* (q, T)) 

Gn((?,r) =n((7°) (e^'^°G(g,T) + G*((?,T)) (33) 

So, for positive 

Gi 2 (q, T) = 2n(g°)Re G{q, T) = ^ Re Gn (g, T) (34a) 

With a small change of notation, G{q,T) = iF{q,T), this becomes 

-fFi 2 {q,T) = 2n{q^) Im. F{q,T) = ^ Im Fu{q, T) (346) 

This relation between the (12) and (11) elements of F{q,T) is reminiscent of the zero-temperature 
unitarity relation, which expresses the imaginary part of the forward elastic scattering amplitude 
(which is the Fourier transform of a matrix element of a time-ordered product of fields) to the sum 
of integrated squared amplitudes needed to calculate reaction rates . So here the imaginary part of 
Fii[q,T) (which is the Fourier transform of a thermal average of a time-ordered product of currents) 
is related to Fi 2 iq,T), which up to a factor of i is the sum (30a) of integrated squared amplitudes 
needed to calculate a thermal reaction rate. 

At zero temperature the amplitude is real for q^ < 4m^, and for q^ > 4m^ twice its imaginary part 
is the discontinuity across the set of branch cuts that run along the positive real axis from q^ = Am? 
to infinity. Such a relationship holds for each Feynman graph that contributes to the amplitude; 
the Feynman-graph relationship is known as the cutting rule, and is named after Cutkosky who hrst 
discovered it^^^]. 

There are also cutting rules at nonzero temperaturel^^. The simplest graph that contributes to F{q, T) 
is _ 

(thermal graph) 


This graph has a branch cut that runs from = 4m^ to oo, like at zero temperature, and another 
from = 0 to —oo. It is real between these cuts. The discontinuities across each of these are 
proportional to the imaginary part of Fii{q,T). However, this fact is not very useful, because for 
higher-order graphs there is no gap between the branch cuts and no region where Fii{q,T) is real: it 
cannot be analytically continued from the upper side of the cuts to the lower side, and so the notion 
of a discontinuity is not useful. The relation (34b) between Fi 2 {q,T), from which the reaction rate 
is calculated, and the imaginary part of Fii{q,T) is still valid, but it is no longer the discontinuity 
across a cut. 

When we calculate Fii{q,T), the internal lines of the graph correspond to thermal propagators Dii{k) 
and Dii[q—k). From the explicit form (21a) of Du, we see that these internal propagators are singular 
at k? = m? and {q — k)'^ = m?. When we perform the necessary integration over the internal loop 
momentum k, these singularities generatet^^^ the branch points at = 0 and 4m^. However, Dyson 
resummation replaces the bare propagators Du with dressed propagators 
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so that the graph becomes 


(thermal graph) 


The resummation removes the singularity of Dii{h) at kP' = m? and replaces it with a singularity 
whose position, according to (25), is given by 

P = m2 + n(A:,r) or + n*(A:, T) (35) 

The consequence is that the branch point of Fn that was at = 4m^ before the resummation has 
now moved off the real axisI^^K But along the real axis the relation (34b) between Fi 2 {q,T) and 
the imaginary part of Fii{q,T) is still valid. 



Infrared divergences 

The infrared divergences of zero-temperature become much worse at finite temperature: the Bose 
distribution diverges at zero energy and causes the usual logarithmic divergences to become power 
divergences. We know that the infrared divergences must cancel if the theory is to make sense, and in 
practice they always do, but there is no general theory yet to show this. To some extent, the situation 
can be rescued by including thermal self-energy insertions in the propagators, so that they acquire a 
mass proportional to the temperature. But, in the case of photons or gluons, not all the degrees of 
freedom have a mass, according to perturbation theory. If the so-called magnetic mass is non zero, it 
is nonperturbative. 

Consider, for example, the effect on the decay rate —> e“’"e“ of the microwave background, which 
is a heat bath consisting only of photons^^^l. Let T be the decay rate in vacuum. The heat bath will 
change it partly because it gives the electrons an additional temperature-dependent mass Sm^ oc e^T^. 
This causes a change dT = Sm‘1 dT / dm^ , which is associated with thermal graphs of the form 


There is also the thermal graph 


One finds that 




^ ^ + J d^k5{k‘^)n{k^) j - ml)5^^\p^ - ml) 

, N 2 

( ^ FP2-P)- (Pi + P2 + A: - P)} 

\Pi-k p2-kj ^ 

where = {\ml — ml). In the integral, the first -function corresponds to the contribution from 
the internal vertices in the thermal graphs being both 1 or both 2 and the second to 12 and 21. Each 
term separately is infrared divergent, like f dk/k"^, but the divergences cancel. 

In fact there is more cancellation than just that of the infrared divergences. For T <C Q one can 
expand (36) in powers of One finds that the first term, of order oemT"^, exactly cancels 

the electron-mass-shift contribution dT, so that the net change in the decay rate is of order uemT^- 
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The lesson is that, in thermal field theory, it is of importance to calculate all terms; gauge theories at 
finite temperature are rife with cancellations. 


Another example where there is an apparent infrared problem is that of the calculation of the equation 
of state for the quark-gluon plasma. For instance, in the purely gluonic case thermal graphs of the 
form 


..'I 




I \ 
I I 
I 


/ 


become more and more divergent as more vertical lines are added. However, the problem goes away if 
one sums over all thermal graphs^^^^. Consider scalar field theory, where the Hamiltonian H contains 
a term f Insert this into Z = tr and differentiate with respect to m^: 

T^logZ = -i I £xmx)f) = {{m?)V (37a) 


Here I have used translation invariance; {{(j){x)Y) is independent of x. With the relation PV = Flog Z 
from (4) and the definition (22c) of £>( 2 ( 3 ;), this gives 


dP 

dm? 



1 

2 


d^q 




(376) 


Insert the form of D'^ 2 {q) obtained from the matrix product D' = MD M in (24) and use the fact 
that the pressure must vanish when the particles are infinitely heavy; 


P = - 



dm'^ 


dS <Qo) y _i_ 

( 27 r)^ — 1 q2 _ ^2 _ 


(37c) 


where H is the self energy defined in (24b) and €{qo) = ±1 according to whether go is positive or 
negative. As any divergence of H now appears in the denominator, the summation has made it 
harmless. One has to worry about the denominator possibly vanishing when g = 0 and m = 0, but 
this will be rendered harmless by the g^ appearing in d'^q = q^dqdVt. The quantitities in (37) are all 
unrenormalised; there is still some work needed, involving the theory of composite operators, to cast 
the formula into one that involves only finite renormalised quantitities!^^]. 


Linear response theory 

Suppose that the thermal equilibrium of a plasma is disturbed by the switching on at t = 0 of an 
external electrostatic potential Ag^j.(x). Then the system’s Hamiltonian acquires an extra term 

H'{t) = e{t) I d^xJ\x)Al^,{x) (38) 

where is the charge density. This will cause to change. As it is a Heisenberg-picture operator, 
its equation of motion is 

^ = i[H + H', jO] (39a) 

where H is the original Hamiltonian. When we take the thermal average of this equation, the con¬ 
tribution from H will disappear because originally there was thermal equilibrium. So the integrated 
change in {J^{x)) at very large time is 

6(J°(x)) = J d'^x'Gn{x — x')A^^^(x') (396) 
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where 

Gr{x - x) = 9{t - t'){[J°{x), J°(x')]) 


Taking the Fourier transform, 


6{J^{k)) = Gn{k)Al,ik) 


(39c) 

(39d) 


We may express the retarded Green’s function Gr in terms of elements of the matrix Green’s function 


Gr 


f/^00 /^OO I /^OO /^00\ 

1^11 “ ^22 + Lt 21 “ ^ 12 ) 


(40) 


Each of the terms here may be calculated from pertubation theory. However, it is simpler to express 
Gr in terms of the function G'^^lk) that appears in the diagonal matrix associated with (see 
(23)): 

Gr = n(A:°)e5^l^“l{(G°° + G°°*) sinh + (G°° - G°°*) sinh (41a) 


On the other hand 


Gii = n(A:°)e^^l''“l{(G°° +G°°*)sinhi/3|A;°| + (G°° - G°°*) cosh i/3fc°} (416) 


So it is sufficient to calculate G^i and change its imaginary part by multiplying it by tanh ^(3k^. 
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